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Abstract 



The paper presents new method for calculating the low-temperature asymp- 
totics of free energy of the 3D Ising model in external magnetic field {H ^ 0) . 
The results obtained are valid in the wide range of temperature and magnetic 
field values fulfilling the condition: [1 — tanh(/i/2)] ~ e, for e C 1, where 
h = f3H, (3 - the inverse temperature and H - external magnetic field. For 
this purpose the method of transfer-matrix, and generalized Jordan- Wigner 
transformations, in the form introduced by the author in H, are applied. 



PACS number(s): 05.50. +q 
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I. FORMULATION OF THE PROBLEM 



As is well known, till now an exact solution for the 2D Ising model in external magnetic 
field (H 0) was not found. In the case of the 3D Ising model there does not exist an 
exact solution forvanishing magnetic field (H = 0), without even mentioning the case with 
magnetic field. Despite great successes in investigations of the Ising model reachedusing 
the renormalization group method ||, and other approximate methods [|3|-[J, the problem 
of calculation various asymptotics for the 2D and 3D Ising models in external magnetic 
field (H ^ 0) is still of great importance. In the paper [|7| we calculated low-temperature 
asymptotics for the 2D Ising model in external magnetic field (H ^ 0), and free energy for 
this model in the limit of asymptotically vanishing magnetic field. In this paper we shortly 
discuss the problem of calculation of the low-temperature asymptotics for free energy in the 
3D Ising model in external magnetic field (H ^ 0), following the approach and the ideas we 
introduced in the paper 0. 

Let us consider a cubic lattice built of N rows, M columns and K planes, to vertices 
of which are assigned the numbers <J nm k from the two-entries set ±1. These quantities we 
will be calling here and everywhere below the Ising "spins." The multiindex (nmk) numbers 
vertices of the lattice, with n numbering rows, m numbering columns, and k numbering 
planes. The Ising model with the nearest neighbors ineteraction in external magnetic field 
is described by the Hamiltonian of the form: 

NMK 

H ^ ] ( s J\@nmk@n+l,rnk J2@nmk@n,m-\-l,k + J3&nmk&nm,k+1 + H&nmk) , (1-1) 

(n,m,fc)=l 

taking into account anisotropy of the interaction between the nearest neighbors (t/1,2,3 > 0), 
and the interaction of the spins cr nmk with external magnetic field H, directed "up" (a nm k = 
+1). The main problem consists of calculation of the statistical sum for the system: 

Uh)= E - E e-^ = 

<thi=±1 (TiVMK=±l 



E ex p 

{"nmfc=±l} 



.nmk 



, (1-2) 



where i^i,2,3 = (3Ji,2,3, h = /3H, (3 = 1/ksT. Typical boundary conditions for the 
variables a nm k are the periodic ones. We follow this standard assumption everywhere below. 
Let us note here that the statistical sum (1.2) is symmetric with respect to the change 
(h^-h). 

In this letter we consider a limited version of the problem. Namely, the problem of 
calculation of the low-temperature asymptotics for free energy in the 3D Ising model in 
external magnetic field. More precisely, given the coupling constants (</i,2,3 = const) and 
external magnetic field (H = const), we consider the region of temperatures satisfying the 
condition: h ~ e" 1 , £ < 1. To be more exact, we introduce a small parameter in the 
following way: 

1 - tanh(/i/2) ~ e, e < 1, (1.3) 
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Then we consider the problem of calculation of free energy per one Ising spin in the ther- 
modynamic limit, with exactness up to quantities of the order ~ e 2 in expansions of the 
operators associated with interaction of spins as well among themselves as with the external 
field, (details of the approximation used will be presented below). In our opinion the prob- 
lem formulated above is of reasonable importance, and, as far as is known to the author, it 
was not investigated in the existing literature. 



II. PARTITION FUNCTION 

Let us consider an auxiliary 4D Ising model in external magnetic field H on simple 4D 
lattice (N x M x K x L). We write the Hamiltonian for the 4D Ising model with the nearest 
neighbor interaction in the form: 

J'2®nmkl@n,m-\-\,kl J?>@nmkl@nm,k+l,l Ji^nmkl^nmk,l+1 ~\~ HtJnmfcij 

n,m,k,l 

(2.1) 

taking into account anisotropy of the interaction between the nearest neighbors (Ji,2,3,4 > 0), 
and interaction of the spins a nm ki with external magnetic field H, directed "up" {pnmki = 
+ 1). Here (2.1) the multiindex (nmkl) numbers the vertices of the 4D lattice, and the 
indices (n, m, k, I) take on values from 1 to (N, M, K, L), respectively. As in the case of the 
3D Ising model, we introduce periodic boundary conditions for the variables a nm ki- Then 
we write the partition function Z A {h) in the form: 



Mh)= E - E <r /m = E exp 



E {Kl a nmklO'n+l,mkl + 



<Tim=±l cr NMKL =±l {c„ m H=±l} Lnmkl 

K2(TnmklO'n ! m+l,kl + K3&nmkl(?nm,k+l,l + K \o 'nmklO 'nmk ,1+1 + hanmkl)} , (2-2) 

where the quantities Ki and h are defined as above (1.2) |||9|]. The expression (2.2) we can 
write, using the well known method of transfer matrix, in the form of a trace from the L-th 
power of the operator T: 

Z 4 (h) = Tr(f) L , f = T^T^T^Tl 12 , (2.3) 

where the operators ^2,3,4,^ are defined by the formulas: 

T\ = exp I K 1 E T nmk T n+l,mk ) > T 2 = exp I K 2 E T nmk T n,m+l,k ) 5 (2-4) 
V nmk / \ nmk / 



T 3 = exp (k s E <mk<mMi) > T ± = ( 2 sinh 2K A ) NMK ' 2 exp (k* 4 E <mk) , 

\ nmk / \ nmk / 



(2.5) 



T h = exp (h E r z nmk J , (2.6) 

\ nmk / 
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and the quantities K4 i K\ are coupled by the following relations: 



tanh(i^ 4 ) = exp(-2Kl), or sinh 2K 4 sinh 2K* A = 1. 



(2.7) 



The Pauli spin matrices r^jf commute for (nmk) 7^ 
matrices satisfy the usual relations 
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n'm'k'), and for given (nmk) these 
It is easy to see that the matrices Ti^^h commute 
among themselves, but do not commute with the matrix T4. In the case in which one of the 
quantities Ki = 0, (i = 1, 2, 3), we get obviously the known expressions describing the 3D 
Ising model on a simple cubic lattice. Namely, the transition to the 3D Ising model with 
respect to the coupling constants K\, K 2 , or K 3 is realized by taking (K x = 0), or (K 2 = 0), 
or (K 3 = 0), and removing summation over n, (N = 1), or over m, (M = 1), or over k, 
(K = 1), respectively. As a result we get the standard expressions || for the 3D Ising 
model in external magnetic field. In the process the operators Tj, (i = 1, 2, 3) in every one 
of the cases are identically equal to the unit operator (T, = 1). A bit different situation 
appears in the case of the transition to the 3D Ising model with respect to the coupling 
constant K±. In this case we take (K 4 = 0, L — 1), i.e. we remove summation over /. In 
consequence we get the following expression for the operator T 4 , (2.5): 



T 4 * = T 4 (fT 4 = 0) 



n< 

nmk 



1 + r, 



nmk 11 



(2.8) 



where we used the relation (2.7). Then, after transition to the limit (K4 = 0, L = 1) 
in (2.3), we can write the following expression for the partition function for the 3D Ising 
model: 



Z 3 (h) 



Tr(TlT]l 2 T z T 2 T x T^) 



(2.9) 



where the matrices Tj are defined as above (2.4 — 6,8). Now we pass to the fermionic 
representation. For this aim one should write the matrices Tj in terms of the Pauli operators 

± IFJl. 



r, 



nmk ' 



T, 



nmk 



it: 



+ 1T V \ 

ik nmk) 1 



(2.10) 



which satisfy anticommutation relations for one vertex, and which commute for different 
vertices. 

As the next step one should pass from the representation by Pauli operators (2.10) to 
the representation by Fermi creation and annihilation operators |T|. In the paper |I[ were 
introduced appropriate transformations (generalized transformations of the Jordan- Wigner 
type), enabling the transition to the fermionic representation: 



T nmk eX P 



N M k-1 



N m-1 



n-1 



X X/ ^sp^spk + X/ a smk a smk | a nmk 



1~nmk eX P 



T nmk ~ eX P 



s=lp=lq=l s=l p=l s=l 

(N M k-1 n-1 M m-1 

EI PspqPspq + EI EI PspkPspk + EI PnpkPnpk 

s=lp=lq=l s=l p=l p=l 

N m-1 K N k-1 n-1 
'~ I EI EI EI 

itpqlspq itmqlsmq + EI Ismklsmk 



iS=l p=l q=l 



=lq=l 



=1 



nmk 



Inmk 
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JV m-1 K n-1 K fc-1 

T nmk = eX P ^ | E E E VspqVspq + E E VsmqVsmq + E VnmqVnmq 
s=l p=l q=l s = l q=l q=l 

n-1 M K M fc-1 m-1 

s=l p=l 5=1 p=l g=l p=l 

n-1 M ftT m-1 K fc-1 

T rimfc = eX P m | E E E @tpq@spq + E E ®npq@npq + E ®nmq@n 
s=l p=l q=l p=l q=l q=l 



Vnmk 



'nmq 



CO 



nmk 



nmk 



(2.11) 



and analogously for the operators r~ mk . In the paper Q we obtained formulas for relations 
between various Fermi operators, and commutation relations for them. Further in this paper 
we will use the fact that the following equality of local occupation numbers is valid: 

T nmk T nmk = a nmk a nmk = finmkfinmk = Inmklnmk = VnmkVnmk = ^nmk^nmk = ^nmk^nmk- (2-12) 

Then, applying the expressions (2.10 — 12) and considerations from the paper 0, we can 
write the partition function (2.9) in the form: 

Z 3 (h) = {2cosh 2 h/2) NMK < 0|T*|0 >= A < 0\U + fi 2 CUD\0 >, U = T l h T 3 T 2 T{T r h , (2.13) 

where A = (2 cosh 2 h/2) NMK and \x = tanh(/i/2), and the operators T 12 ,3, T^ r and C, D are 
of the form: 

N,M,K 

Ti = exp Kt ( a nmk ~ «nmfc)(«n+l,mfc + «n+l,mfc) 

ra,m,fc=l 

N,M,K 

T 2 = exp K 2 E (Pnmk - Pnmk)(Pn,m+l,k + Pn,m+l,k) 
n,m,fc=l 

N,M,K 

T 3 = exp K 3 (®nmk - 9nmk)(0nm,k+l + ®nm,k+l) J ( 2 -14) 

n,m,fc=l 



and 

T h = ex P j/^ 
T l h = exp \ fx 2 



N-n 



Al-m 



EE 



nmk n+s 



,mk ^ y ^ y "nmfc^n', 



nmk s=l 

K-k 



nn'mk t=l 

M-m 



nn'mm'k 1=1 
N-n 



^ ] ^ ^ Qnm,k+l@nmk ^ ] ^ ] &n,m+t,k^nmk' ^ ] ^ ] @n-\-s,mk@nm' k' 

.nmk 1=1 nmkk' t=l nmm'kk' s=l 



(2.15) 



C — E ^nmfc, -P — E 



nmk 



Here and below Snm... means summation over the complete set of indices (n = 1, ..JV; m 
1, ...M; etc.). It is obvious that the operator G: 



G 



S E a nmk a nmki 
nmk 



(2.16) 
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where S is the operator of the total number of particles, commutes with the operator T*, 
(2.13). Therefore, we can divide all states of the operator T* into states with even (Ag = +1) 
or odd number of particles (A^, = —1) with respect to the operator G, (2.16). The form 
of the operators T 12 ,3 does not change during the course, only the boundary conditions for 
the operators (aw*:,---) do. In the case of even states (A^ = +1) antiperiodic boundary 
conditions , and in the case of odd states periodic ones, are chosen [[/]]. 
The next step is transition to the momentum representation: 

exp(i-K/4:) ^ -i(nq+mp+kv) t+ f)+ . n + Q+ _> 

/ j Smw; h'nmk Iqpvi v nmk So 



(iV M ' Ky^fo ^qpui t-'nmk ' 'Iqpvi "mnk ' ^>qpvi 

and introduction for fixed (qpv) corresponding bases for 77- and (- Fermi creation and 
annihilation operators in the representation in terms of occupation numbers in the finite- 
dimensional Fock space of dimension 2 8 = 256). Then, after a series of transformations and 
calculations we arrive at the following formula for the partition function (2.13): 

M=4 II W) ( II ^(y)) <0\T*(h)T 2 Tl(h)\0>, (2.17) 
where the operators T*(h), T 2 , T 3 *(/i) are of the form 

-^1 (^) ~~ eX P ^ 1 (?) {£,-q-p-v£qpu £-q-pv£qp-v £-qp-v£q-pv £-qpu£q-p-h 

0<q,p,v<n 

T 2 = exp I 2K 2 l C0S P(vtpuVq P » + ■■■) + sinp(r)± q _ p _ v r}+ v + ... + riq P vV-q- P -» + ••• 

I 0<q,p,u<w 



T£(h) = exp 



y "j B^iy^i^qppQ—q—p—u ~\~ Q—qpvCq—p—v ~t~ Cq—puQ—qp — v + C—q—pvCqp—l 
0<q,p,u<ir 



(2- 



and Ai(q, h), ... are defined by the expressions: 

Ai(q, h) = cosh 2/^ — sinh2K 1 cosg + a(h, q) sinh2ii' 1 sing, 
A$(v, h) = cosh 2K3 — sinh 2K% cos v + a(h, v) sinh 2K 3 sin u, 
a(h, q) [cosh 2K\ + sinh 2K\ cos q] + sinh 2Ki sin q 



B 1 (q,h) 
B 3 (u, h) 



Ai(q,h) 

a(h, v) [cosh 2K% + sinh 2K3 cos v] + sinh 2K% sin v 
A 3 (u,h) ' 



a(h, q) = tanh 2 (V2) 1 + C ° Sg , a(h, u) = tanh 2 (/i/2) 1 + C ° SZ/ . (2.19) 

sin q sin v 

In the formula for Z^ D (h) the sign (+) means that we consider the case of even states 
(Ag, = +1) with respect to the operator G, (2.16). It is obvious that for h = we arrive at the 
3D Ising model in vanishing magnetic field. Then, for K\ = (or K 2 = 0, or K 3 = 0) the 
expression (2.17) for the statistical sum describes the 2D Ising model in external magnetic 
field m. 
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III. SOLUTION 



Let us consider calculation of free energy per one Ising spin in external magnetoc field 
in the approximation described shortly in the introduction. For this aim let us consider the 
operators T*(h) and T£(h) in the "coordinate" representation: 



T*(K) = exp 
T*(h) = exp 



N-n 

a (. S ) a nmk a n+s,mk 

nmk s=l 

K-k 

E E <i)e 

nmk 1=1 



(3.1) 



where the "weights" a(s) and c(l) are defined by the formulas: 



as 



c(i) 



^ £ 2fi 1 (g)sin( S g) = 2:r + tanh 2 at 7 i- 

1 



K 



2 J B 3 (z/) sin(/z/) = z* 3 l + tanh 2 h 



z- 



0<i/<7T 



We introduced renormalized quantities (if* 3 , h\. 



defined as follows : 



1,2,3, 



/ = 1,2,3, 



(3.2) 



sinh 27^*3 = /3 13 [sinh 2^1,3(1 - tanh 2 (/z/2)], 



cosh(2ir* 3 ) = /5i i3 [cosh2K 1>3 + tanh 2 (/i/2) sinh2.Fr : 



1.3jj 



/3 li3 = [1 + 2tanh 2 (V2) sinh2K 1 , 3 e 2 ^- 3 ]- 1 / 2 , tanh 2 a? 3 = tanh 2 (a/2) - /?1 ' 3 ex ^ 2Kl ^ 



(3.3) 



These formulas are valid for 3 > 0). As in the case of the 2D Ising model [J7|, LTJ], also in 
this case one can introduce a diagrammatic representation for the vacuum matrix element 
S =< 0|T 3 *(a)T 2 T 1 *(a)|0 >. Computation of the vacuum matrix element S, which enters the 
formula (2.17) for Z 3D (h) in general case, where the "weights" (3.2) are arbitrary is, at least 
at present, impossible. Nevertheless, there exists a special case in which we can calculate 
the quantity S in the 3D case. Namely, this is the case where the "weights" (3.2) are 
independent of / and s. In this case one should, as in the 2D case []I2" |, put the parameters 
Ki t s equal zero (i£i )3 = 0) in the formula (2.13), and then express the operators Tjf in terms 
of the Fermi /3-operators (2.11) of creation and annihilation, with the goal to calculate S. 
After transition to the momentum representation, one should calculate the vacuum matrix 
element S*(y x , y 3 , z 2 ): 



S*( yi ,y 3 ,z 2 ) =< 0\T\y 3 )T 2 T r ( yi )\0 >, y 1>3 



tanh 2 h 



1,3, 



where z 2 = tanhi^ 2 becomes trivial. (Here we introduced the following change of notation: 
h/2 — > hi - for the operator T[, and h/2 — > h 3 - for the operator T l h ). We can write the 
result for S*(yi,y 3 , z 2 ) in the following form: 

S*( yi ,y 3 ,z 2 ) = (2cosh 2 fG)^ JT [(1 - 2z 2 cos p + z%) (I - cosp) + 2z 2 ( yi + y 3 ) sin 2 p+ 

0<qpu<ir 

1 4 

2/12/3(1 + 2z 2 cosp + z 2 )(l + cos p) . (3.4) 
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This result can be used further to calculate free energy in the approximation discussed above 
(1.3). For this aim let us note that the conditions [tanh 2 h\ 3 /(l — z\ 3 ) 2 ] — > 1 are equivalent, 
accordingly to (3.3), to the conditions (exp(— 2i^i j3 )(l — tanh 2 /i/2) — > 0). It follows from 
this equation that for fixed (Ji,3 = const, H = const) these conditions are satisfied in the 
region of temperatures T, in which (h/2) ~ e" 1 , e <C 1. In this case we can use the 
result (3.4). Namely, let us consider the formulas (2.19) for -61,3, written in terms of the 
renormalized parameters (h\ 3 , K13)' 

tanh 2 ht smgH _|_ 2 Z * s j n q( u \ 
1 - 2^ 3 cosg(z/) + 2* )3 

where 2* 3 = tanh if* 3 . Next, since the following equalities are satisfied: 

4,3 21,3(1 -tanh 2 /i/2) 



1 + ^/ 1 + 2z 1)3 tanh 2 h/2 + z\ 



3 



then, if we introduce a small parameter [1 — tanh(/t/2)] ~ e, (e <C 1), and expand B\ r 
into a series in powers of e (z* 3 ~ e), we obtain 



(tanh 2 /z* 3 + 22* 3 )sing(z/) 

-di,3 = ; r- ~ e 

1 — cosg(z/j 
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This formula gives the following expressions for the " weights" a(s) and c(l), (3.2) in this 
approximation: 

a(s) = tanh 2 /i* + 2z{, c(l) = tanh 2 h* 3 + 2z* 3 , (3.6) 

with exactness of the order of smallness ~ e 2 . As a result in this approximation the "weights" 
a(s), c(l) do not depend on (s,l). Finally, if we substitute to the expression (3.4) for 
S* (yi, V3, Zi ) the parameters y\ — > a(s) and y 3 — ► c(/), (3.6), we arrive at the following 
formula for free energy on one Ising spin F^o{h) in the thermodynamic limit: 

- (3F 3D (h) x ln(2 3/2 cosh K{ cosh K 2 cosh K* 3 cosh 2 h/2) + — /* In [(1 - 2z 2 cosp + z\) x 

27T JO L 

(1 - cosp) + 2z 2 (tanh 2 h\ + tanh 2 h* 3 + 2z{ + 2z 3 ) sin 2 p + (tanh 2 h\ + 2^*)(tanh 2 h* 3 + 

2z 3 )(l + 2zi cosp + -2 2 )(1 + cosp) dp, (3 

where (3 = l/k B T, and z 2 = tanhif 2 , and h\ 3 and K^ 3 are coupled with h and i^i )3 by the 
relations (3.3). One can show that, as it was done for the ID and 2D Ising models [0,0, in 
the case of the states odd (A^ = —1) with respect to the operator G, (2.16), the formula for 
F 3D {h) is described in the thermodynamic limit by (3.7). Let us note that the asymptotics 
(3.7) obtained above can be applied also in the case of rather strong magnetic fields (if), as 
far as it satisfies the condition (1 — tanh h) ~ e, e -C 1, (T — const). 
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IV. FINAL REMARKS 

The result derived above(3.7) can be applied to the analysis of equilibrium thermody- 
namics of the three dimensional Ising magnetic, lattice gas, and also three dimensional 
models of binary alloys [L4j] in the region of temperatures and magnetic fields (1.3) de- 
rived above. Such analysis, as well as construction of appropriate phase diagramms for the 
models mentioned above is, in our opinion, of great interest. They deserve presentation in 
a separate publication. Therefore we deliberately do not compare here our result (3.7) with 
the existing papers devoted to this problem. The other important feature of the presented 
method is the possibility of deriving expressions for the free energy of the 3D Ising model in 
the limiting case of the magnetic field tending to zero (H — > 0, N, M, K — > oo) if we know 
exact solution for the 3D Ising model in the absence of external magnetic field (H = 0). 
This possibility results from equations (3.2 — 3.3) describing renormalised interaction con- 
stants If* 3 , and corresponds, as was presented in the paper J/J, to the results obtained by 
C.N.Yang |0| for the 2D Ising model. 
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